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Chapter 14 Gauss’ theorem 


We now present the third great theorem of integral vector calculus. It is interesting that 
Green’s theorem is again the basic starting point. In Chapter 13 we saw how Green’s theorem 
directly translates to the case of surfaces in R? and produces Stokes’ theorem. Now we are 
going to see how a reinterpretation of Green’s theorem leads to Gauss’ theorem for R?, and 
then we shall learn from that how to use the proof of Green’s theorem to extend it to R”; the 
result is called Gauss’ theorem for R”. 















































A. Green’s theorem reinterpreted 


We begin with the situation obtained in Section 12C for a region R in R?. With the 
positive orientation for bdR, we have 
| ty 


j | Shavay 7 
ye sededy = ~ f ode. 


bdR 














We now consider a vector field F = (F), F2) on R, and we obtain from Green’s theorem 


bdR 


(Notice if we were thinking in terms of Stokes’ theorem, we would put the minus sign on the 
left side instead of the right.) 
We now work on the line integral, first writing it symbolically in the form 


F e (dy, —dr). 
bdR 
We now express what this actually means. Of course, bdR may come in several disjoint pieces, 
so we focus on just one such piece, say the closed curve y = 7(t) = (y(t), 72(t)), a< t <b. 


Then the part of the line integral 
corresponding to this piece is actually 


[ FoW) «hb, -ah(0ya 
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This in turn equals 





[ FOU) we ln’ @)llae. 


The unit vector in the dot product has a beautiful geometric interpretation. Namely, it is 
orthogonal to the unit tangent vector 


(4), 2) 
Iv @l 


and it is oriented 90° clockwise from the tangent vector. (You can use the matrix J of 


Section 8D to see this, as 
! eer, 
1-2) 
V2 V1 


You can also think in complex arithmetic: —i(y;, + 175) = 74 — iy{, and multiplication by —i 
rotates 90° clockwise.) 

Because the tangent direction 7/(t) has the region R placed 90° counterclockwise from it, 
we conclude that the normal vector 


(y(t), ~7i(¢)) 
yl 








points away from R: 


n~ 


N 
Let us name this unit normal vector N. Thus at each point (x,y) € bdR, the vector N (x, y) 
at that point is determined by these requirements: 





e N has unit norm, 
e Nis orthogonal to bdR, 
oN points outward from R. 


Green’s theorem has now been rewritten in the form 


/| (f+ a) ae dy = [ Pe Ravon, 
bdR 
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This result is precisely what is called Gauss’ theorem in R?. The integrand in the integral 
over R is a special function associated with a vector field in R?, and goes by the name the 
divergence of F’: 














Again we can use the symbolic “del” vector 
0 600 
: - (s =) 
Ox’ Oy 


divF = VeF. 


to write 


Thus Gauss’ theorem asserts that 


/ | VeFdvoly = | F e Ndvol,. 
R 
R 


Notice that this final result has the interesting feature that the orientation of the coordinate 
system has completely disappeared. 





B. Gauss’ theorem for R” 


It is very easy now to imagine what the correct extension to R” should be. This is very 
different from our encounter with Stokes’ theorem, for which the final result seems to be very 
much tied to R%. 

Suppose D C R” is an open set with a “reasonable” boundary bdD. We don’t want to 
explain this further, except to say that at each point of bdD there is a normal vector N such 
that 











e N has unit norm, 
e Nis orthogonal to bdD, 
e N points outward from D. 


We are willing to allow corners and edges to some extent, but nothing too “wild” is to be 
permitted. 
Suppose F is a vector field of class C! defined at least on clD, the closure of D. We define 
the divergence of F by the formula 
OF, OF, 


ivF = Fs —-4..-4+——, 
div Ve aa, t + ar 
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This is again the formal dot product 


O O 
VeFf = (goers ger) # Fie Fa) 


GAUSS’ THEOREM. Under the above hypotheses, 





[ve Pavot, = [ Fe Nave. 
7 bdD 











We shall spend the remainder of this section discussing examples of the use of this theorem, 
and shall give the proof in the next section. 

First, apply the theorem to the very particular vector field F(x) = x, so that Ve F =n. 
Then we obtain immediately 


nvol,(D) = / re Ndvolp_1. 
bdD 
In particular, if D is the ball B(0,r) of radius r, we know vol,(D) = a,r". And bdD = the 


sphere S(0,7), for which N= x/r. Thus xe N= ||2||?/r =r, so we obtain immediately from 
Gauss’ theorem the equation 


Voln—1(9(0,a)) = nayr”’. 


This formula was found with more effort in Section 11C. 
For another example we examine an ellipsoid. Let A be a symmetric positive definite n x n 
matrix, and D the “solid” ellipsoid 











D = {x ER" | Arex < 1}. 





We know from Problem 10-42 that 


An 


vol.) = vere 





We also have the gradient 
VArexr = 2Az, 


so that in Gauss’ theorem 
A Arex 1 
Te — = 


Axl] Aa 
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Thus we obtain 





/ dvoln-1 Ay 
| Az Vdet A 
bdD 





In particular, if A is the diagonal matrix with entries aj”,...,a7?, we obtain 
dvoln—1 
= Ana,...An, 

xy wR, 

cee ap tt ad 

A 4.-+ B=1 : . 
a2 a2, 
1 n 


a result obtained in Problem 11—28 with much more effort required. 


PROBLEM 14-1. Prove that if f is a real-valued function and F is a vector field, 
then the product rule is valid: 


Ve(fF) = fVeF+VfeF. 











PROBLEM 14-2. Prove these special relations for R°: 





0, Le. 
0, 


divergence of a curl 


Ve(V x F) 


and 


Ve(FxG) = GeVxF-FeVxG. 


Perhaps the most important use of the Gauss theorem is that it affords us a geometric 
interpretation of divergence. This interpretation is in the same spirit as our previous dis- 
cussions of gradient, in Section 2F, “Geometric significance of the gradient,” and of curl, in 
Section 13G, “What is curl?” 

We approach this by first understanding the geometry of the integral over bdD in the right 
side of Gauss’ theorem. The terminology that is used derives from applications to fluid flow 
or electric fields. If F’ is a vector field in R” and M is a hypermanifold for which a preferred 
“side” has been chosen by naming a continuous unit normal vector N, then F'e N is the 
component of F' orthogonal to M, with a definite sign affixed. Then the integral 














| Fe Ndvoln—1 
M 
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in a sense measures the net amount of the vector field that “crosses” M. This is usually called 
the net flux of F' across M. 
Then the right side of Gauss’ theorem, the integral 


/ FeN dvoln—1, 
bdD 


represents the net flux of the vector field in the direction outward from D across its boundary. 
Thus Gauss’ theorem becomes in words, 


i VeFdvol, = flux of F' out from D. 
D 





Now let 29 € R” be a fixed point, and choose D to be the ball B(x, €) with center xo and 
radius ¢€. If € is small, we have a good approximation 


Z V e Fdvol, ~ V e F(x) X voln(B(Zo, €)). 
B(ao,€) 


Thus we obtain the formula 
1 


ee a F e Ndvoln-1, 
Ve F(xo) <0 Voln(B(o, €)) hon — 


which we express in words as 
VeF(x) = the flux of F away from 29, per unit volume. 


Hence the word “divergence.” 
Once again, we have an amazing situation that although V e F’ was originally defined in a 
way that crucially required a coordinate system, the fact is that the quotient 


1 | 7 
es es Fe Ndvoln_1 
voln(B(X0,€)) J s(2o,€) 


manifestly is independent of any coordinate system; therefore, its limit as « — 0, namely the 
divergence of F at 29, is also completely determined geometrically. 


PROTOEXAMPLE. This example in case n = 3 is the force field of a unit electric charge 
placed at the origin, or the gravitational field of a unit mass at the origin. In the case n = 3 
it’s called the “inverse square law.” For general n, we have to within a constant factor 


x 


F(a). = —— 


|x" 
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Then we calculate from the product rule of Problem 14-1, 


Ver = ae +V((lzx||") #2 





n me 
= jee * (ll 2r)er 


n = 
jer el mall? 


0. 





Of course we notice the crucial role here of the power n. It comes from the fact that on R” 
we have Vex =n. Thus on R” our vector field is 


ra) — ill 


flere? 








so it’s an “inverse (n — 1) power law.” 


Now we compute an easy case, the outward flux of F’ across the sphere S(0,a). As the 
outward unit normal is N = «/a, we find 


i. Fe Ndvoln_1 = | = * dvoly_1 
x/a S(0,a) Sr ee 


= vol,_1(.9(0, 1)) 


= Nn. 


NOTATION. In the present investigations the number na, is going to be appearing quite 
often. Let us denote it as 


Wn = NAy 
= vol,i($(0,1)): 





That is, w, is the measure of the unit sphere in R”. Thus 
Wg = 27, 
W3 = An 5 


Wy = Qn, 
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and in general 
Qn”/? 
= =e . 11-15). 
w T(n/2) (see p ) 
Now suppose D is a region in R” such that the origin is exterior to it; that is, 0 € clD. 
Since F’ is well defined in clD and has zero divergence, Gauss’ theorem implies immediately 
that 





| F e Ndvoln_1 = 0. 
bdD 


O. 





On the other hand, suppose the origin is in the interior of D: 0 € D. Then we can prove 
that 


/ Fe Ndvoly_1 = Ws. 
bdD 


This is seen very easily by using what is often called a “safety ball.” That is a ball B(0,€) 
whose radius is so small that clB(0,¢€) C D. 





Then D — B(0,e) is “safe” in that the origin is not in it, so we have from the preceding case 


/ Fe Ndvoly_1 = 0. 
bd(D—B(0,e)) 


But bd(D — B(0,€)) is composed of the two disjoint portions bdD and the sphere S(0,€). 
Thus 


/ Fe Ndvoly_1 +f Fe Ndvoly_1 = (). 
bdD S(0,€) 


The latter integrand has N = —2/e, as N is oriented outward relative to D — B(0,e). Thus 
the latter integral is actually equal to —w,,. This proves the desired equation. 
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Another viewpoint of the above analysis is to move the center of the field from the origin 
to an arbitrary point 2 by translation. Then we obtain 


/ Og ee 
bap || — Zol|” 0 if a ¢clD. 


The case zp € bdD does not fall in either category. 


PROBLEM 14-3. One might guess that for 79 € bdD 


= ~ | 
/ cece e Ndvol,_1 = <Wp- 
pap ||z — Zol| 2 


Show that this is indeed the case if n = 2 and 3 and D is a “ball” (disk if n = 2). 


PROBLEM 14-4. Show that the result of the preceding problem is valid for all n > 2. 


NAMES. In Russian texts Gauss’ theorem is called Ostrogradski’s theorem. It’s also frequently 
called the divergence theorem. 


C. The proof 


In one sense the proof we give is a generalization of the proof of Green’s theorem as given 
in Chapter 12. The basic idea of integrating first in the direction of the partial derivative in 
the integrand is still the key ingredient of the proof. That is, we want to prove that for each 


fixed index 7 
OF; 


p Ox; 





dvol, = / FN e €;dvol,_1, 
bdD 


and so we approach this by performing the x;-integration on the left side before the other 
integrations. We then sum for 7 = 1,...,n, to obtain the theorem. In our analysis we may as 
well replace the component F; with just a real-valued function f of class C'. Thus we need 
to prove 





| of avo, = f Nidvoln_1, (x) 
p O02; bdD 
where N; = Ne e; is the 7** component of the outer unit normal vector N. 

In another sense the proof we present here is much more efficient than the one we gave in 
Chapter 12, as we employ a beautiful technical device to get around the troublesome details 
about patching results together (see Section 12C). This technique is called a partition of unity. 
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Here is a brief description of it. Suppose € > 0 is given. Then there is a set of C1 functions 
U1,-..-,Wwn defined on R” such that 





N 
Souk = 1 on clD, 
k=1 


for each k, w, = 0 outside some ball of radius e. 


This is easier to prove than you might imagine. One first chooses nonnegative w,’s with 
the second property, such that their sum is positive on clD, and then each is normalized by 
dividing by their sum. Incidentally, instead of C! we could use C? or C® etc., or even C®. 

Using such a partition of unity, it suffices to prove (*) for each product vf and then add 
the results to achieve the theorem. Another way of expressing this is to say that we may 
assume ab initio that f itself is zero except on some small ball. 

There are now two separate cases to handle. The first is the case in which the ball where 
f is nonzero is contained in (the interior of) the open set D. Then the right side of (*) is zero 
since its integrand fN; equals zero on bdD. But also the left side is zero, as we can write 











Of Of 
LL. = dvoly 
[ oe ors vO 
Pub i ( il = te, dvOln1 
R-1 —oo Ox; 
and the inner integral 
de, = 0 
=e Ox; an 


by the fundamental theorem of calculus. 

That was easy. The second case is not much harder but is very much more interesting. In 
this case f is nonzero only in a small ball, but this ball intersects bdD. Here we explain the 
philosophy behind our partition of unity. We choose the balls so small that when one of them 
intersects bdD the open set D near that ball has a particular sort of description. Namely, 
it can be represented as “one side” of a graph of a C! function. For instance, consider the 
following typical situation: 
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The portion of D in the illustrated ball is described by an inequality of the form 
Be OO ican ee reyes cla 

(If the ball were on the “right side” of the picture, the defining inequality would appear as 
ORG gs 2A): 


We content ourselves with the first case in the proof.) 

We pause to streamline some of the notation. First, since f is zero outside the small ball, 
we may as well assume f is defined for all x; satisfying the given inequality. 

Second, we may as well arrange coordinates so that 7 = n. This merely has the effect of 
greatly shortening the notation. 

Third, we denote by x’ = (x,...,%pn—1) so that f is defined for all points x satisfying 


te, > ye), 2 eR, 











where R C R"* is some small ball. 7 
We now have formulas for N and dvol,_;. Namely, since N is the outward unit normal, 


N _ (Oy /On,, ee ,00/O%n-1, =1) 


7 JVel? +1 | 








and 
dvol,_1 = J||Voll? +1 dey we. dLn-1 
— Jer ai ae 
Thus 
Nes. = ae . ae. 
—dz' i t=, 


These formulas require us to finish with two different cases. 








z= n| In this easier case we have 


O Fay Fubini ! 
oe [ a ) Be Dds 

















FTC [rl da! 
In=y(z’) 
- [tee 
R 
— fNndvoly_1, 


bdD 
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and (x) is proved. 














i#n| Here we have to face the problem that doing the x; integration first may not be so 


feasible. We need to do the x, integration first, so for each fixed x’ € R we shall change 
variables by setting 


Then 0 < t < co and we have 


oi dvol, "2" | | le uae 
p Ox; RJ g(a’) OXi 
rJo O 


Poe i. of —(x', p(x’) + t)dtda’ 
R 


x (0,00) ) OX; 


chain rule 0 ! nN | Op Of ! 
mn i Ee (F(a', o(a") + 8)) — FE | dae’ 








(x’) + t)dtda’ 








We now have two integrations to perform. We handle the first by doing the x; integration 
first. The result is zero, as for each fixed t we essentially have 
= : ; FTC 
— (f(a', p(x") +#)) de; "EP 0. 





In the second integral we return from ¢ to x,. Thus we have 








of Op Of / 
n= ; t)dt 
[peeve i © a Baa p(x’) + t)dtdx 


= By lemipeenrs 
[-#f- on ot ee 


Lpn=0o 








rro f Ov, , 





da’ 


R In=(x') 


0 
= f SEF@ oar) ae! 


= f Nidvol,_ 1- 
bdD 








This finishes the proof of Gauss’ theorem. 
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D. Gravity 


The applications of Gauss’ theorem are quite diverse and important. Especially in the 
rudiments of electricity and magnetism and in gravitation we find the theorem playing an 
important role. We content ourselves here with some observations on gravitation. 











Isaac Newton’s theory includes as an axiom that two point masses in R® attract one another 
with an inverse square force. We of course are determined to generalize this to R", where it 
presumably makes no physical sense. So suppose we have a mass m, located at 2 € R” and 
another one mz located at y € IR”. Then the force at y due to the mass at x is equal to the 
vector 














Xx — 
Gimymz Pedal 

I|z — || 

where the number G > 0 is called the “gravitational constant.” We shall assume from now on 
that G = 1. 


One of the first things that can be proved has to do with a mass distributed uniformly 
over a sphere. We can set this situation up in the following way. Let us suppose the sphere 
is centered at the origin, so that it equals S(0,a), where a is its radius. Let us suppose the 
mass spread over this sphere is 7. We say the mass density is the constant 


M M 


VOoln_1(S(0,a)) wat)" 





Now suppose that we have a unit point mass located at y. Then according to Newton the 
force at y due to the mass spread over the sphere is given by the integral 





ae M 
| = = az dvoln_1(2). 
S(0,a) Iz — y||" wna 


What we are going to do now is calculate this integral. 


In doing the calculation we may choose coordinates to suit ourselves, so let us place the 
positive n™ coordinate direction pointing toward y. In other words, we assume that y = ré,, 
where of course r = ||y||. 
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By symmetry we notice that the above integral, which is a vector, must have zero components 
in the directions orthogonal to é,. Thus we may as well just compute the n coordinate of 


the integral, that is 
M Le 


wpan! S(0,a) lz — ré,||” 


THE PHYSICAL CASE n = 3. It is quite impressive that the easiest dimension for performing 
the above integral is n = 3. We now carry this out. First, we use the familiar result 


dvol,_1(2). 


Iz — yl? = Jal? — 22 ey + |lyll? 


to rewrite the denominator, obtaining 


M 3-7 
Ana? Js(o,q) (a? — 2rx3 + r?)3/? 








dvolo(z). 


Next use the usual spherical coordinates with 73 = acos y, to obtain 


M 9 [ acosp—Tr ae aa 
-Qr a” sin 
Ara? (a? — 2ar cos p + r?)3/2 oe 








0 
=f acosp—Tr 
2 Ja (a? — 2ar cos yp + r2)3/2 


sin ydy. 





PROBLEM 14-5. Physics texts know what to do here. Make the substitution 
u=a- 2ar cos (p +r? 


to change the dummy ¢ into u, and thus show that the above integral equals 


—M/r? ie a SG, 
—M /2a? i oS 
0 if O<r<a. 


The physical interpretation of this result is well known and extremely important. Namely, 
at points outside the closed ball B(0, a) the gravitational attraction of the uniform mass spread 
over the sphere is the same as if all the mass were concentrated at the center of the sphere. 
At points in the interior B(0, a) of the ball, the gravitational attraction is zero. And at points 
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on the sphere $(0,a) itself, the gravitational attraction is the same as if all the mass were 
concentrated at the center of the sphere, divided by 2. 


PROBLEM 14-6. (not recommended) Repeat all of the above for the case n = 2. 
The integration is much more difficult; a table of integrals might be helpful. 


At this point you should be wondering what any of this has to do with Gauss’ theorem. 
We now explain this. The direct calculation of the integral on the preceding page is quite a 
formidable task. The case n = 3 is not so bad, but the case n = 2 is extremely difficult, and 
one can imagine that the situation for general n could be even worse. Now we are going to 
show how a combination of some beautiful Euclidean geometry and Gauss’ theorem enable us 
to calculate the integral very easily. 


INVERSION. The geometry we have just referred to is called inversion in a sphere. (See 
pp. 6-34 and 7-23.) Suppose that S(0,qa) is a sphere in R” and y € R”. As you will see, we 
must assume y # 0. We then say that the inversion of y is the point y* which is uniquely 
defined by requiring that it lie on the ray from 0 through y and that ||y||||y*|| = @?. In other 
words, 








2 
a 
= a ae Uh 
Ile? 
Notice that 


and that 


* 


y = y—y e S(0, a). 


PROBLEM 14-7. Verify the following picture: we assume ||y|| > @ and that a line 
has been drawn from y tangent to the sphere at the point z. The picture is then in the 
unique two-dimensional plane determined by 0, y, and z. 
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PROBLEM 14-8. As in the preceding problem, verify this similar two-dimensional 
depiction of the relation between y and y’*. 


(Om 


Now we notice a very interesting phenomenon. If x € S(0,a), then 


lla —y* |? = lial? — 22 © y* + |y* |)? 
; a’y at 


— 27 e —— + —— 
llyll?— llvll 


= 7p (Ilyll? -— 22 ey +a’) 








a 3 (Ilyll? — 22 y + [lzl|’) 














=i ally — all’. 


That is, 





lc —y"|| = FT lx — y|| for all x € S(0,a). 
y 


This equation is quite interesting, in that it shows that the ratio ||x — y*||/||a — y|| is constant 
in its dependence on x € S(0,a). 

We now return to the calculation of the integral at the bottom of p. 14-12. We of course 
may assume y # 0, since when y = 0 the integral equals zero thanks to the oddness of 
the integrand. By symmetry it is clear that the vector represented by this integral has zero 
components in directions orthogonal to y; in other words, 


M ZY 


—___. ———— dvoln_i1(z) = —cy, 
wna"! Js(o,a) lz — yll” 
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where c is a scalar. In other words, the gravitational attraction at y points toward the origin. 
As we just need to calculate c, let us form the dot product of each side of the equation with 
y, denoting r = ||y||. We obtain 


M pey—r 
aay cue Sec dvol,_1(z) = —er?. 
wna” Jgo,a) lle — yll” 


As we are wanting to apply the divergence theorem, we would like the integrand to be presented 
as the dot product of a vector field with N, which is equal to x/a. Therefore we write the 
numerator 


2 
crey-r = rey—— rex 
a 


I| 
TN 
< 
| 
8/3 

iw) i) 
8 
WF 
e 
8 


Aha! Notice how naturally the inversion y* has appeared! We thus have found that 


M = ~ 
c= —— —— e Ndvol,-1(2). 
wna” Js(o,a) lle — yll 


Now we utilize the relationship between ||x — y|| and ||2 — y*|| to find that 


M ae 
c= — | aa a e Ndvolp_i(2). 
wnt Js(o,a) lle — y*l| 


We find that we are in the situation of the protoexample of Section B. There are three 
cases. 





\|y|| >a] In this case ||y*|| < a, so we obtain from p. 14-8 that 











M 
C= Hs 
rn 
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\|y|| <a} In this case ||y*|| > a, so we obtain instead 











c = 0. 





|y|| =a] In this case ||y*|| = a, so Problem 14-4 gives 











M 


2a” 


In summary, the gravitational attraction of a uniform mass distribution over a sphere in 
IR” behaves in the exterior of the sphere as if all the mass were concentrated at the center, 
and is zero in the interior of the sphere. At points on the sphere itself, it is as if all the mass 


were concentrated at the center, divided by 2. 





E. Other differentiation formulas 


PROBLEM 14-9. In analogy with the formula for curlF’ given on p. 13-22, suppose 
that {%1,.--, Qn} is an orthonormal basis for R”, and represent points in coordinates 


relative to this basis as 














i t1P1 oe te ba Pie: 


Then explain why 
“0 
ivf = —(Fe &j). 
div d. a, | e (i) 


PROBLEM 14-10. — Explain why the formula just presented can also be written in 
terms of directional derivatives as 


divF = \° D(F eG; Gi). 
a=] 


Then show that this formula remains true for a general orthogonal basis. 





DEFINITION. Given a vector A € R", the first-order differential operator Ae V is defined 


by 
(AeV)(f) =AeVf, 


where f is a real-valued differentiable function. 
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REMARK. Since both gradient and dot product have intrinsic meaning, so does this oper- 
ator. In fact, it is clear that in terms of directional derivatives 


(AeV) f(x) = Df (a; A). 





PROBLEM 14-11. Show that for vector fields F and G on R® we have 











Vx(F x G)=(VeG)F-(VeF)G+(GeV)F—-(FeV)G. 














PROBLEM 14-12. Show that for a vector field F on R?, 
Vx(Vx F)=-WF+V(VeF). 
In this formula the Laplacian of F' is the vector field 


3° OF 


VF = . 
= Ox? 





F. The vector potential 


This section corresponds to Section 12F, where we considered the problem of finding a 
(real valued) potential for a vector field with zero curl. That is, given a vector field F on R? 
satisfying V x F’ = 0, we sought a potential function f for which F = Vf. We based our 
thinking on the differentiation formula 





VxV7=0 


(see Section 13G). 
The analogous formula we now investigate is the one which states 


Ve(VxG)=0 


for any C? vector field G on R? (see Problem 14-2). We are thus naturally led to the following 
question: if F is a vector field on R® for which V e F = 0, does there exist a vector field G 
satisfying 








FHV XG? 
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The answer is yes, provided that F is a vector field defined on all of R*. The proof amounts 
to a rather logical procedure of integrating the above differential equation for G. Before doing 
this, note that there is a strong lack of uniqueness. For if F' = V x G, and also F' = V x Go, 
then 0 = V x (G, — G2). Thus from Section 12F we conclude that G; — Gp = Vf for some 
function f. Conversely, since V x (Vf) = 0, we can always add to G any gradient field without 
changing the validity of F = V x G. 


DEFINITION. A vector field G satisfying the equation 
F=VxG 
is called a vector potential for the field FP’. 


THEOREM. Suppose F is a vector field of class C' defined on all of R°, and suppose 
Ver =0. Then there exists a vector potential for F'; this vector potential is unique up to the 
addition of a gradient field. 





PROOF. The differential equations which G must satisfy are 


OG3 OG, 

i esc Yay 
Ox» 0x3 = 
OG, _ 3G) _ 
0x3 Ox iii, 
dG, OGy . 
Or, — Ax = Fs. 


This task seems daunting at first glance. However, we may as well assume G3 = 0, for we can 
find a function f whose gradient equals 


Vi = (?, ce G3) 


by the simple device of writing 


Then G — Vf has zero third component. 
With G3 = 0, the equations now become 


SUES i 
Ars — fl; 
OG, 
ebay 
Ors 2) 
OE ae 
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Now choose G2 to be any function satisfying the first of these equations. It remains to solve 
for G. Integrate the second equation to find G; to within an additive function f (2x1, x2): 


Gi = Gi + f(21, 22). 
Then we need to find f so that the third equation is satisfied. That is, 
Of __p G9 . OG 


Oty Oty Oy 
There is no problem integrating this equation, provided that its right side is independent of 
x3. We check this easily: 
O 0G? AG. 
Ors 4 * Bay" a 
OF3 7 0G? ra OG 
0x3 §=Ox30%2 = =OX302X, 
_ OF 0G? OG 
7 Ax; 7 02x20X3 7 0x102x3 
OF; OF, OF, 

Or, «OX ORY 

= (by hypothesis. 











QED 


PROBLEM 14-13. Show that the above procedure produces for a constant vector 
field F the vector potential G = (#3F — x2F3, —23F},0). 


PROBLEM 14-14. If F is constant we earlier calculated that SF x £ is a suitable 
vector potential (see Problem 13-10). Find a function f so that the solution of the 
preceding problem satisfies 


1 


PROBLEM 14-15. Using the notation of Problem 13-11, suppose F(Z) = AZ is a 
linear vector field. Prove that V e F = traceA. Supposing traceA = 0, find a vector 
potential G for F’. 
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PROBLEM 14-16. Suppose f and g are real valued functions of class C? defined on 
an open set D C R®. Define the vector field 
F=Vf x Vg. 


Prove that V e F = 0, and find a vector potential for F’ which is defined on D. 
(HINT: compute V x (fVqg).) 


PROBLEM 14-17. Show that the vector field 


F= 0 
(ae e+ y? 


has zero divergence on R*— (z-axis). Show that it has a vector potential of the form 
f(a, y)k, and calculate a suitable f explicitly. 





PROBLEM 14—-18*. Show that the vector field 


x y 
Fe= 0 
(start ) 


has zero divergence on R®— (z-axis), and find a vector potential for F' which is defined 
on all of R?— (z-axis). 





PROBLEM 14-19. (The analog of Section 12G.) Suppose that F' is a C' vector field 
defined on a star shaped open set D C R?. Suppose that F has zero divergence. Supposing 
that D is star shaped with respect to the origin, define the vector field 


1 
| Pia) < tedk. 
0 


Prove that 
VxXGeHF. 
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PROBLEM 14-20. Use the preceding problem to give another solution of Problem 
14-15, and show that the corresponding vector potential is 


1 
3(Ad) x a. 


PROBLEM 14-21. Consider the gravitational vector field 


(x,y, z) 
(a + y? + #?)3/2 





on R?— origin. It has zero divergence. Use the scheme of the general theorem of this 
section to find a vector potential. Show that you may reach an answer in the form 
Zz 
G => GP —2, 0). 


(a? +9? )a/a? + ¥? + 2? 





PROBLEM 14-22. The solution which you obtained for the preceding problem may 
seem unsatisfactory, as it is not only badly behaved at the origin (necessarily), but also 
on the entire z axis. Prove that in fact there is no vector potential for the gravitational 
field F defined on all of R°— origin. 

(HINT: supposing such a potential G would exist, compute 


/ F e Ndarea 


unit sphere 





two different ways.) 


